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A METHOD FOR COMPUTATION OF VIBRATION MODES 
AND FREQUENCIES OF ORTHOTROPIC THIN SHELLS OF REVOLUTION 
HAVING GENERAL MERIDIONAL CURVATURE 

By Howard M. Adelman, Donnell S. Catherines, 
and William C. Walton, Jr. 

Langley Research Center 

SUMMARY 

This report describes a procedure for computing the vibration modes and frequen- 
cies of thin shells of revolution having general meridional curvature and orthotropic elas- 
tic properties. The procedure is based on the finite- element method in which the direct- 
stiffness approach is used. A geometrically exact finite element is employed. A com- 
puter program based on this procedure has been written and details of the program are 
described. The geometric characteristics of the shell are used as inputs to the program 
in the form of functions of the meridional coordinate. The stiffness and mass matrices 
are computed by numerical integration by use of the trapezoidal rule. 

The computer program is applied to several shell configurations including two 
cylinders, two conical frustums, shells of both positive and negative Gaussian curvature, 
and an annular plate. Frequencies are correlated with frequencies from previous inves- 
tigations for these shells. The agreement between results of the present analysis and 
results from the previous investigations is generally excellent. 

INTRODUCTION 

A problem of current interest to structural analysts in the aerospace field is that of 
determining the dynamic behavior of structures in which some of the components are thin 
shells of revolution. Understanding the modes of vibration of the individual shell com- 
ponents can be of fundamental importance in connection with this problem. Consequently, 
much effort has gone into developing techniques to determine natural frequencies and mode 
shapes of the commonly encountered shells of revolution. Few closed-form solutions are 
known and, therefore, most of the developments have been in the area of approximation 
methods. Among the methods that have been tried are Rayleigh- Ritz methods (refs. 1 
and 2), Stodola-type iteration methods (ref. 3), finite-difference solutions (ref. 4), finite- 
element methods (refs. 5 to 9), and methods in which the shell boundary-value problem is 
reduced to an initial-value problem involving first-order differential equations which 


are numerically integrated (ref. 10). The finite-difference and numerical- integration 
methods involve a trial-and-error search for the natural frequencies that will make a 
certain determinant vanish. These "search methods" are relatively slow, and analysts 
using them have been known to overlook modes, as noted in reference 3. Stodola-type 
methods also lose numerical significance in the calculation of higher modes as noted in 
reference 3. 

The authors, in the course of developing practical procedures to analyze the forced 
response of structures incorporating shells of revolution, required a method for computing 
mode shapes and frequencies of shells of revolution having general meridional curvature 
and orthotropic elastic properties. These mode shapes and frequencies would be used in 
analyses of structures involving such shells where in the analysis the deformation of each 
shell is represented by superposition of a number of mode shapes of the shell. Early 
experience indicated that selection of representative modes for a shell would require 
examination of a great number of its modes some high in the frequency spectrum. It was 
therefore necessary that: 

(1) The method should give capability for quick calculation of a large number of 
modes and frequencies 

(2) The mode shapes and frequencies quite high in the frequency spectrum should be 
accurately predicted 

(3) The analyst should be protected from overlooking modes in computation 

In view of these objectives, search methods and Stodola methods were considered 
unsatisfactory for the reasons of their inadequacy to meet these requirements. Both the 
Rayleigh-Ritz and finite-element approaches seemed to offer better chances for success 
in meeting the objectives. 

From the viewpoint of the analyst, the outstanding advantage of finite-element and 
Rayleigh-Ritz approaches is that they lead to a symmetric eigenvalue problem which is 
amenable to fast and accurate solution on a digital computer. In the methods available 
for solving symmetric eigenvalue problems, all the modes are computed simultaneously, 
and thus any danger of overlooking modes in computation is avoided. 

A finite- element approach was selected in preference to a Rayleigh-Ritz approach 
for the following reasons: 

(1) The computing details of the Rayleigh-Ritz methods reported in references 1 
and 2 resulted in use of a large number of terms. These methods lead to relatively poorly 
conditioned eigenvalue formulations requiring that a large number of significant figures 
be carried in the calculations in order to retain significance in the results. Preliminary 
trials with circular plates indicated that the finite- element approach leads to very well- 
conditioned eigenvalue problems. 
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(2) It was believed that the finite- element method would converge more easily than 
the Rayleigh-Ritz method when local high stress gradients, such as occur for some edge 
conditions (ref. 11), are present. 

The element most popularly employed in the finite -element analysis of shells of 
revolution has been the conical element (for example, refs. 8 and 9). This element can 
exactly fit cylinders and conical frustums. However, for shells having a curved meridian, 
use of this element leads to only an approximation of the shell by a series of joined conical 
frustums. Thus, the curved meridian is approximated by a series of straight lines. Con- 
sequently, an analysis of a shell with a curved meridian based on conical elements may 
give inaccurate frequencies and stresses (refs. 5 and 6). 

Analysts have been aware that the use of an element which coincides with the shape 
of the shell being analyzed would probably improve the accuracy in computed results 
(ref. 5). The main impediment to the use of such a geometrically exact element has been 
a reluctance on the part of analysts to give up a certain computational convenience asso- 
ciated with the conical element. This convenience is that since the shape is fixed, quadra- 
tures required to compute the stiffness and mass matrices of the element are performed 
only once, and the same matrices are used in every analysis. With a geometrically exact 
element, the shape of the element depends on the shape of the portion of the shell which the 
element represents, with the result that the integration has to be an inseparable part of 
each analysis. It has been recognized that a natural and probably feasible approach to 
making the quadratures part of the analysis is to use numerical integration (ref. 5). How- 
ever, the objection has remained that for each element the radius of the shell and the two 
radii of curvature must be specified as functions of position along the meridian of the 
element. 

In spite of this objection, the decision was made to develop a computer program to 
meet the previously stated objectives based on a geometrically exact element. It was 
believed that the necessity for description of the geometry of an element in terms of func- 
tions rather than of numerical parameters would present no difficulties in practice if, as 
is nearly always the case, the geometry of the entire shell could be described by functions 
located in a subroutine which could be readily changed. Development has progressed to 
the point where the program has been applied to a variety of shells of revolution of practi- 
cal interest. Detailed correlations have been made between frequencies from the program 
and frequencies calculated for these shells by other investigators. A cursory correlation 
of mode shapes including stresses has been made for some of these shells but as yet is 
inconclusive because of present unavailability of sufficient modal data from the methods of 
the previous investigators. These correlations are not presented in this report. 

The main purposes of this report are as follows: 

(1) To describe the analysis underlying the computer program 
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(2) To describe the computer program 

(3) To present the frequency correlations 

SYMBOLS 


a 0 k’ a 2 ^ coefficients in polynomial displacement function for normal dis- 
placement w 

matrix which transforms displacements and rotations at ends of an element 
to coefficients of polynomial displacement functions (see eq. (13) and table I) 

b Q ^jbj j^bg ^b^ k coefficients in polynomial displacement function for meridional 

displacement u 

B matrix defined in equation (42) 


c 0,k’ c l,k’ c 2,k’ c 3,k 


coefficients in polynomial displacement function for circumferential 
displacement v 


Ck matrix whose elements are coefficients in an expression for strain energy of 

a shell element in terms of coefficients of polynomial displacement functions 
(see eq. (20)) 

Cn,Ci 2 ,C 22 membrane stiffness constants 

Cgg in-plane shear stiffness 

D diagonal matrix whose elements are eigenvalues of mass matrix (see eq. (39)) 

Dh,Di 2 ,D 22 flexural stiffness constants 

Dgg torsional stiffness 

mi ddle-surface strains (see eqs. (la) to (lc)) 

Ek kinetic energy of kth element 

E kinetic energy of shell; also Young's modulus 
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frequency 


f 

matrix whose elements are coefficients in an expression for kinetic energy of 
a shell element in terms of coefficients of polynomial displacement functions 
(see eq. (30)) 

h shell thickness 

I identity matrix 

i = fi 

K number of elements used to represent a shell 

K 11’ K 12>^22» K 66 stiffness constants representing interaction between in-plane and 

out- of -plane strains 

L meridional length of a shell 

m meridional wave number for a freely supported cylinder 

Mk element mass matrix 

M shell mass matrix 

n circumferential wave number 

N order of stiffness and mass matrices after edge constraints have been applied 

P matrix whose elements are coefficients in an expression for kinetic energy 

of a shell element in terms of displacements u, v, and w (see eq. (28)) 

q index representing an integration station 

Q total number of integration intervals 

r radius of a shell measured in plane normal to shell axis 

Rf,R2 principal radii of curvature of shell 
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R 

matrix whose elements are coefficients in an expression for strain energy of 
shell element in terms of actual variables in strain energy (see eq. (16)) 

> 


s 

meridional coordinate 

s k 

element stiffness matrix 

s 

shell stiffness matrix 

s 0 

meridional distance from origin of s to reference edge of a shell 

s k 

meridional distance from reference edge of shell to center of kth element 

s kq 

meridional distance from reference edge of shell to qth integration station 
of kth element 

s ki 

meridional distance from reference edge of a shell to ith location on kth 
element at which mode shape is evaluated 

t 

time 

T k 

inverse of matrix 

u 

meridional component of middle-surface displacement 

U 

matrix whose columns are eigenvectors of mass matrix 

V 

circumferential component of middle-surface displacement 

v k 

strain energy of kth element 

V 

strain energy of shell 

w 

normal component of middle- surface displacement 

X 

meridional coordinate measured within a single element (see eq. (5)) 

x kq 

meridional distance from center of kth element to qth integration station 
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X matrix which describes assumed form of variables appearing in strain energy 

y column matrix containing unknown displacements and rotations 

Y matrix which describes the assumed form of displacements u, v, and w 

Z matrix defined in equation (43) 

/3 rotation of shell generator relative to unstrained direction (see eq. (10)) 

6 modal column (see eq. (44)) 

y column matrix whose elements are coefficients of assumed- displacement 

polynomials (see eq. (18)) 

meridional length of kth element 

6 circumferential coordinate 

Kl,K2> K i2 changes in curvatures (see eqs. (Id) to (If)) 

XpXg,. . .,X^ eigenvalues of mass matrix 

I column matrix whose elements are displacements and rotations at ends of an 

element (see eq. (22)) 

p mass density 

u) circular frequency 

fi nondimensional frequency 

p Poisson's ratio 

Primes denote differentiation with respect to s or x; superscript T denotes 
transpose of a matrix. 

Special notations used in machine plots of figures 16, 17, 20, and 21: 

N circumferential wave number 

S/L nondimensional meridional distance 
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u,v,w 


middle- surface displacements in the meridional, circumferential, and normal 
directions, respectively 


UMAX,VMAX,WMAX maximum values of U, V, and W, respectively 

= (Circular frequency)^, sec”^ 

DEVELOPMENT OF THE STIFFNESS MATRIX FOR A 
GEOMETRICALLY EXACT ELEMENT 

Strain Energy in Terms of Displacements 

For purposes of the following analysis, reference is made to figure 1. In this fig- 
ure, u, v, and w represent displacements in the meridional, circumferential, and 
normal directions, respectively, Rj and R 2 are the two principal radii of curvature 
of the shell, and r is the radius of the shell measured in a plane normal to the shell 
axis. All three radii are regarded as functions of the meridional coordinate s, measured 
along the shell from a reference edge. 

According to Novozhilov (ref. 12), the six strain- displacement relations which 
describe the local state of strain for a thin shell of revolution are as follows: 

Membrane strain in meridional direction: 



(la) 

Membrane strain in circumferential direction: 

1 9v 1 , w 

2 r 90 r R 2 

(lb) 

In-plane shear strain: 

e 12 = ~ I 7 + v * ” - r ’ v 

12 r 90 r 

(1c) 

Change of curvature in meridional direction: 

K, = -w" + u’ — Ri ’u 

1 R 1 R X 2 1 

(Id) 

Change of curvature in circumferential direction: 

_ 1 9^w 1 9v r'w’ 1 , 

2 r 2 qq2 rR 2 90 r + rRj r U 

(le) 

Twist of the middle surface: 

„ 1 9w' , 1 . 9w 1 9u 1 , 1 , 

“12--F1FV '^ + rR 1 90 + R 2 V "rR 2 rV 

(If) 
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For a shell which, in general, is composed of orthotropic layers, the strain energy 
is given in reference 13 (p. 45) as follows: 

V = 2 j J( C ll e l 2 + 2C 12 e l e 2 + C 22 e 2 2 + C 66 e 12 2 ) r de ds 


|II( D 


11 K 1 2 + 2D 12 k 1 k 2 + D 22 k 2 2 + D 66 /C 12 2 ) r d0 ds 


+ JJ(^ll e l K l + K 12( e l /C 2 + e 2 K l) 

where in equation (2) the integrations are taken 
definitions hold: 


+ K 22 e 2*2 + K 66 e 12 /f 12j r d0 ds 

over the shell surface and the following 


(1) C„, c 12’ ^22 are membrane stiffnesses 

(2) Djj, Dj2, D 22 are ffexural stiffnesses 

(3) Cgg is the in-plane shear stiffness 

(4) Dgg is the torsional shear stiffness 

(5) Kn, K 12 , K 22 , Kgg are stiffnesses due to the interaction between in-plane 
strains and changes in curvature 

All of these stiffnesses are, in general, functions of the meridional coordinate s. Ref- 
erence 13 contains an excellent discussion of the derivation of the above stiffnesses for 
shells having various numbers of layers and composed of materials having various types 
of elastic properties. 

The work in the present study is based on Novozhilov's strain-displacement rela- 
tions (eqs. (la) to (If)), the energy expression of equation (2), and the definitions of the 
stiffnesses in reference 13 with the following single exception. The strain (called 

r in refs. 12 and 13) is defined by Ambartsumyan (p. 25) to be double the value of this 
strain as defined by Novozhilov. Since the authors prefer to use Novozhilov's definition 
of /Cj 2 > the value of Dgg used herein is four times the value of Dgg given in refer- 
ence 13 and the value of Kgg is twice the value of Kgg given in reference 13. 

For a shell of revolution vibrating in a natural mode with circular frequency u>, 
the three displacements u, v, and w can be expressed as follows: 


u(s, 0,t) = u(s)cos n0e* w *" 
v(s,0,t) = v(s)sin n0e* wt > 
w(s,0,t) = w(s)cos n0e iwt 


( 3 ) 
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The displacements from equations (3) are substituted into the strain-displacement 
relations of equations (1). Substitution of the resulting strains into the strain- energy 
expression of equation (2) and integration with respect to 6 yields the strain energy in 
terms of displacements. The amplitude of the strain energy for n # 0 is as follows: 



For n = 0 the strain energy as given by expression (4) should be doubled. The suc- 
ceeding developments are carried out on the assumption that n 0 with the under- 
standing that for n = 0 appropriate expressions should be doubled. 

Representation of the Shell by Geometrically 
Exact Finite Elements 

The present analytical method follows the main steps of conventional finite-element 
analysis. It is noted, however, that each element coincides exactly with a slice of the 
actual shell. Hence, the elements are spoken of as "geometrically exact elements." 

A typical idealization of a shell of revolution is shown in figure 2. Counting ele- 
ments from the reference edge, the following definitions are made: 


K 


total number of elements 


length of kth element, measured along meridian curve of shell 


coordinate inside kth element, measured along meridian from center of kth 
interval so that following relationship holds: 

6 k< <fk 
2 = X = 2 


( 5 ) 
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s k distance along meridian from reference edge of shell to center of the 

kth element 

From the foregoing definitions for x and s^, it follows that 

s = s k + x (6) 

A numbering system has been adopted in which quantites such as displacement, 

e k e k 

derivatives of displacements, and rotations at s = s k - — and s = s k + -^ are indi- 
cated by subscripts k and k+1, respectively. Thus, for example, w k is the normal 

e k e k 

displacement at s = s k - y and u k+1 is the meridional displacement at s = s k + 

Also, it is necessary to have a notation for the radius of curvature Rj at the locations 
e k 

s = s k T — . The symbols k and Rj k+ j represent the respective values. 


Assumed Displacement Field for Element 

As an approximation, the displacements u, v, and w are assumed to have the 
following polynomial forms over the kth element: 

w = a 0,k + a l,k x + a 2,k x2 + a 3,k x ^ 


u = b 0,k + b l,k x + b 2,k x2 + b 3,k x3 / 


v = c 0,k + c l,k x + c 2,k x2 + c 3,k xl 


(7) 




where the a's, b's, and c's are undetermined coefficients. From equation (7) it fol- 
lows that: 

r 

l 0,k 


r ^ 

w 

W T 


W T 

u 

U T 

V 


-m 


*l,k 

l 2,k 

*3,k 

3 0,k 

5 l,k 

>2,k 

3,k 

: 0,k 

: l,k 

'2,k 


3,k 


b 


( 8 ) 
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where 


[x> 


1 x x2 X* 3 0 0 0 

0 1 2x 3x 2 0 0 0 

0 0 2 6x 0 0 0 

0 0 0 0 1 x 

0 0 0 0 0 1 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 


x 2 x 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 


2x 3x 2 0 0 

0 1 x 

0 0 1 2x 3x 


(9) 


Relationship Between Undetermined Coefficients and Displacements 
and Rotations at Ends of Element 

The rotation of the meridian curve relative to the unstrained direction is defined 
as /3 and is given by 


It follows that 


and 


u 




' 3 k =w i-sr 


l,k 


u 


%+l " w k+l R 


k+1 


l,k+l 


( 10 ) 


( 11 ) 


( 12 ) 


Inserting x = - y and x = — into the appropriate locations in equation (8) 
results in the following relationship: 


w k 


a 0,k 

u k 


a l,k 

v k 


a 2,k 

£k 


a 3,k 

u k 


b 0,k 

v k 

w k+l 

= [ A k]< 

b l,kl 

b 2,k f 

u k+l 


b 3,k 

v k + l 


c 0,k 

0k+l 


c i,k 

u k + l 


c 2,k 

_ v k+l_ 




( 13 ) 
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where the elements of matrix jA^j 
the following relationship results: 


are given in table I. When equation (13) is inverted, 


a 0,k 


fVl 

a l,k 


u k 

a 2,k 


v k 

a 3,k 


£k 

b 0,k 


u k 

J b l,k 

| b 2,k 


M 

w k + i? 

b 3,k 


u k + i 

c 0,k 


v k+l 

C l,k 


^+1 

c 2,k 


u k+l 

c 3 k 


w 


(14) 


where 


l-l 


N ■ H' 

The elements of the inverse matrix are given in table II. 

Formulation of Element Stiffness Matrix 
From equation (4) the strain energy of an element may be written as follows: 

V 2 


(15) 


V = l 
V k 2 


r 
w 


w f 

w" 


|w, w', w", u, T > T JIXK U 

u’ 


e k/ 2 


dx 


(16) 


v’ 

V. J 


where [r] is a 7 x 7 symmetric matrix, the elements of which are known functions of the 
meridional coordinate x. The elements of jjtj are listed in the appendix. Using equa- 
tion (8) in equation (16) permits the strain energy to be written in terms of the undeter- 
mined polynomial coefficients as follows: 




V 2 


( 17 ) 
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where 


or 



( 18 ) 


v k = 



(19) 


where 



Finally, use of the transformation expressed by equation (14) gives the strain energy as 


V k = |0) T [Tk] T [Ck][T k ]{{) (2D 

where 


r 


w. 


"\ 


u k 
v k 
P k 
u,’ 


© = 


w 


u 


k+1 

k+1 

v k + l 

%+l 


u k+l 

v k+l 
v _/ 


( 22 ) 
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Inspection of equation (21) identifies the shell element stiffness matrix jSj^j as 

M = C T k] T [ C k][ T k] 


( 23 ) 


DEVELOPMENT OF CONSISTENT MASS MATRIX FOR 
THE GEOMETRICALLY EXACT ELEMENT 


If rotary inertia is neglected and if the shell is assumed to be vibrating in a natural 
mode, the amplitude of the kinetic energy E^ for n 4 0 is 

E k = 2 J ^ / /3h ( u2 + y2 + w2 ) r d* (24) 

_e k / 2 

where ph is the mass per unit area of the shell. The quantity ph is a known function 
of s and therefore of x. As with the strain energy, the expression for kinetic energy 
must be doubled for n = 0. 


Based on the assumed displacements of equation (7), the following relation may be 
written: 



where 


1 

X 

X 2 

X 2 

0 

0 

0 

0 

0 

0 

0 

0 



M- 

0 

0 

0 

0 

1 

X 

X 2 

X 2 

0 

0 

0 

0 

(26) 



0 

0 

0 

0 

0 

0 

0 

0 

1 

X 

x2 

X 2 



Equation (24) can be rewritten in the form 


where 




phr 

0 

0 


0 

phr 

0 


0 

0 

phr 


(27) 


(28) 
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Using equation (25) with equation (27) yields 


where 

Fk '*I £ £ k/ lH T [ p I Y ] dx 

_e k/ 

In view of equation (14), equation (29) may be written as follows: 
Therefore, the element mass matrix may be identified as ^M^j, where 

M ■ M T W[Tk] 


DERIVATION OF MODAL EQUATIONS 


(29) 


(30) 


(31) 


(32) 


In view of the numbering convention adopted for the elements, the second edge of the 
kth element coincides with the first edge of the k+lst element. In this analysis, the fol- 
lowing conditions of compatability are assumed to hold at each such juncture: 


r 

w k + i 


r 

w k + i 


u k+l 


u k+i 


v k+l 

\ = ) 

v k+i 

\ 

%tl 

\ 

%+i 


u k+l 


u k + i 


v k+l 

V. ^ 

kth 

element 

v k+l 

V. J 

k+lst 

element 


for all k < K 


(33) 


Of the six equalities in equation (33), the first four are standard. The last two, however, 
are valid only for shells having continuous distributions of stiffness. 


The total strain energy V and the kinetic energy E may be expressed as follows: 



(34) 
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where V k and E k are given in equations (21) and (31), respectively. If the summa- 
tions of equations (34) are carried out and use is made of equation (33), the strain and 
kinetic energies may be written as follows: 

> (3 

E - 1 w, ® T [“K ?} J 

where 


S stiffness matrix, which is a symmetric positive semidefinite or positive 

definite matrix of order 6(K + 1) 

M mass or inertia matrix, which is a symmetric positive definite matrix of 

order 6(K + 1) 

y a vector containing all of the unknown displacements and rotations 


In the present method, the matrices S and M are constructed by the well-known 
procedure of superimposing element matrices, illustrated in figure 3. As the figure 
shows, the superposition consists of placing appropriate shell element matrices in the 
matrices S and M so that the matrix elements in the lower right 6X6 block of the 
kth matrix add to the corresponding matrix elements in the upper left 6x6 block of the 
k+lst matrix. 


The modal equations for shells with no edge constraints may be derived by mini- 
mizing the quantity E - V with respect to each of the variables in the vector y. This 
minimization is equivalent to the following set of equations, 6(K + 1) in number: 


8(E - V) _ 


9w, : 


= 0 


8 (E - V) = Q 


3u 


9(E - V) _ Q 
3v k 




9(E - V) _ Q 

9 ^k 

= 0 ) (k = 1, 2, . . ., K + 1) (36) 

k 


a(E - V) „ 

" , 


Equations (36) can be expressed in the form 




(37) 
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Rigid edge constraints are incorporated by deleting from the stiffness and mass 
matrices of equation (37) those rows and columns which correspond to displacements and 
rotations that must vanish to satisfy the constraints. The form and character of equa- 
tion (37) are not affected by the deletion of rows and columns from the matrices S and 
M, except that S may become positive definite instead of positive semidefinite. Of 
course, the order of matrices S and M is reduced. 

Equation (37) determines 6(K + 1) natural frequencies and modal columns. The 
computation of these frequencies and modes is a standard operation. However, for com- 
pleteness, some of the reductions involved in the solution are given in the section which 
follows. 

The modal columns consist of values of displacements and rotations at each of the 
K + 1 junctures on the shell. For many purposes a more detailed mode shape is 
required. From equation (7) with x replaced by s - s^, the following equations can be 
written for the mode shape within the kth element: 




w = a 0,k + a l,k( s - s k) + a 2,k( S - S kf + a 3,k( S 

- s *) 3 



u = b 0,k + b l,k( s " s k) + b 2,k( s " s k) 2 + b 3,k( s 

- s *f 

> 

(38) 

v = c 0,k + c l,k( S ' S k) + c 2,k( s “ S k) 2 + c 3,k( s 

- < 



The coefficients aQ ^ through Cg ^ are computed by using equation (14). 



COMPUTATIONAL METHOD AND COMPUTER PROGRAM 


In this section, the details of the computing method are given. Each step is dis- 
cussed and the flow of information is shown in block diagrams. The steps are intended to 
be specific enough to allow a digital computer program to be written by using them as a 
guide. 

The basic organization of the computer program used in the present analysis is 
shown in block diagram 1. 

The input which describes the geometry of the shell and its physical properties is 
first read, and stiffness and mass matrices referred to the coefficients of the displace- 
ment polynomials are then computed. These matrices are transformed so that they are 
referred to the displacements and rotations at the ends of the element. When a stiffness 
and a mass matrix have thus been computed for each element, the matrices are super- 
imposed to form the shell stiffness and mass matrices. If the shell has rigid edge 
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constraints, then appropriate rows and columns are deleted from the stiffness and mass 
matrices to satisfy these constraints. The stiffness and mass matrices then constitute 
the ingredients to an eigenvalue problem which is solved for frequencies and modal 
columns. From the modal columns, the coefficients of the displacement polynomials are 
computed. The detailed mode shapes are then evaluated. 



Block diagram 1. 
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A more detailed block diagram (see block diagram 2) is shown. 


INPUT! 


COMPUTE s k , k=l,2, • • *K 


x 


IPRINT INPUT, s k , AND EDGE CONSTRAINT 


x 


j CYCL E ON CIRCU MFE RENTI AL WAVE NUMBER (n) | 


| COMPUTE SIZE OF STIFFNE SS AND MASS MATRICES | 


| SET STIFFNESS AND MASS MATRICES TO ZERO I 

X 


CYCLE ON ELEMENT NUMBER (k) | 


| COMPUTE TRANSFORMATION MATRIX T,^ 
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Uj CYCLE ON INTEGRATION STATION (^)~| 

[COMPUTE COORDINATES OF INTEGRATION STATI0N~7s kq ~AND x Rq ) I 


X 


COMPUTE MATRICES X(x kq ), Y(x kq ), R(s kq ), AND P(s kq ) 

X 


ICOMPUTE MATRICES SUMC kq AND SUMF kq 



COMPUTE ELEMENT STIFFNESS MATRIX S k AND ELEMENT MASS MATRIX M k 


X 


| ADD MATRIX S k IN PROPER POSITION IN STIFFNESS MATR IX S | 

is > . i — — 


ADD MATRIX M k IN PROPER POSITION IN MASS MATRIX M 



DELETE APPROPRIATE ROWS AND COLUMNS FROM STIFFNESS 
AND MASS MATRICES TO SATISFY EDGE CONSTRAINT 

I 


COMPUTE FREQUENCIES AND MODAL COLUMNS 

i - j 1 

1 


OUTPUT FREQUENCIES AND MODAL COLUMNS.' 



Block diagram 2. 


A detailed discussion of the computing in each block of block diagram 2 is as 
follows: 

BLOCK 1: 


The following input quantities are required: 

K number of elements used to represent shell 

length of kth element, where k = 1, 2, . . K 
Sq meridional distance from origin of s to reference edge of shell 

The following input functions are required: 
l/R^l/l^ reciprocals of principal radii of curvature 
r shell radius measured in a plane normal to the shell axis 

r’,Rj' meridional rates of change of r and Rj 


Cli,Ci2,C22 > C 66 

d 11 > d 12 > d 22> d 66 

k 11 ’ K 12 > k 22> k 66 


stiffnesses 


jjohj mass per unit area 

The following control numbers are required: 


NBEG 


initial value of n 


NLAST final value of n 


Q 


number of integration intervals to be used within each, element 


IPRINT if IPRINT = 0, intermediate matrices are not printed 
if IPRINT * 0, intermediate matrices are printed 

NMODE number of mode shapes to be computed 
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IINT number of locations along each element at which mode shape is to be evaluated 

ICASE edge constraint code (see table III) 


BLOCK 2 : 

S l = S 0 + f e l 

k-1 

s k = s 0 + | e k + 2 6i (k = 2, . . ., K) 

i=l 


BLOCK 4 : 

n = NBEG, NBEG + 1, . . NLAST 
BLOCK 5 : 

KN = 6(K + 1) 

BLOCK 7 : 

k = 1, 2, . . K 
BLOCK 8 : 

The elements of [T^J, a 12 x 12 matrix, are given in table II. 
BLOCK 9 : 

q = 1, 2, . . Q + 1 
BLOCK 10: 


Each element is divided into Q equal intervals for the numerical integration. 
There are then Q + 1 integration stations. The values of x and s at the qth integra- 
tion station of the kth element are, respectively, defined as 


and 


x, = e 


kq k\ q 


q - 1 1 


s, = s, + e 


kq " °k T k\ — Q 


q - 1 1 


(q = i, 2, . 

• Q+ i) 

(q = i, 2, . 

• •> Q + l) 


BLOCK 11 : 

The elements of £xj, [jpj, and [1 r] are given in equations (9), (26), (28), 

and (Al) to (A28), respectively. The matrices Ix^x^Nj and |Y/x kc Jj are found by 
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substituting x kq in equations (9) and (26), respectively. The matrices |p^s k( ^ anc * 
^( s kq) are f° un< ^ t>y substituting s^ into equations (28) and (Al) to (A28), respectively. 


BLOCK 12 : 

fm = ’[ X ( X k<lj T [ R ( S kq)][ X ! x kqj] 

[sUMC k i] = £[c kl ] 

j3UMC kq J . jsUMC k q .jj * [C kq ] (q . 2, 3, . . „Q) 

[suMC k>Q+ ^] = [suMC kQ ] + ±[c k>Q+1 ] 

fkq]-[Y(vf[ P N)][ Y N)] 

[suMF kl ] = £[f^ 

[sUMF kq ] = [sUMF k>q _ l] + [F kq ] (q = 2. 3 ® 

[sUMF kQ+ ^] = [sUMF kQ j + |[F kjQ+1 ] 


BLOCK 13 : 

[c K ] = [sUMC kQ+1 ] 

[F k ]= [SUMF,^ 

[s t ]=[T k ] T [C k ][T k ] 

[Mk] = [T k ] T [F k ][T k ] 

BLOCKS 14 and 15 : 

The manner in which the matrices (either stiffness or mass) are placed in the over- 
all matrix is illustrated in figure 3. 
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BLOCK 16: 


Commonly encountered edge constraints along with the appropriate rows and 
columns to be deleted from S and M are given in table III. 

BLOCK 17: 


Compute by the threshold Jacobi method (ref. 14, p. 397) a modal matrix U and a 
set of eigenvalues Xj, . . ., X N for the matrix M. Then 


0 


x 2 


U T MU = D = 


(39) 


V N 


and 


U T U = I 


where 


N order of matrix M 

I identity matrix of order N 


Since M is positive definite, all diagonal elements of D are positive. 
Compute: 



D 1 / 2 = 


lp2 


|PS 



(40) 


(41) 


B = D -1//2 U t SUD" 1//2 


(42) 
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Compute Z, a modal matrix of B, by the threshold Jacobi method. Then 



0 




( 43 ) 


0 



Z T Z = I 

2 

The values of u> are the squares of the circular frequencies. 
Compute: 


6 = UD" 1//2 Z 


(44) 


The columns of 6 are the modal columns. After the computation of the modal columns, 
insert zeros in the locations which correspond to rows and columns deleted from S 
and M. 

BLOCK 18 : 

The output consists of the following: 

(1) Lists of circular frequency squared o> , circular frequency c o, and 

frequency f. 

(2) For each mode, an array of displacements and rotations as follows: 


Note that zeros again appear in the locations corresponding to deleted rows and columns 
in S and M. 

BLOCK 19 : 

i = 1, 2, . . ., NMODE 
BLOCK 20 : 

k = 1, 2, . . ., K 
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BLOCK 21: 



where the elements of T k are given in table II. 
BLOCK 22: 


For the purpose of computing the detailed mode shapes, each element is divided 
into IINT intervals. The number of locations at which the mode shape is to be evaluated 
is IINT + 1. The value of s at the ith location of the kth interval is defined as s^ and 
is given by 


s ki~ s k + e k 


i - 1 
TINT 


(i = 1, . . IINT +1; k = 1, . . ., K) 


The mode shape over the portion of the shell represented by the kth element is then com- 
puted from: 



The entire mode shape is then constructed by placing the portions end to end. 
BLOCK 23 : 

The mode shape consists of arrays of u, v, and w. 
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APPLICATIONS AND DISCUSSION 


Description of Shells Analyzed 

In order to ascertain the generality and efficacy of the present method for computing 
natural frequencies, a number of applications were made. The following configurations, 
shown in figure 4, were treated: 

(1) An isotropic cylinder with freely supported edges investigated by Arnold and 
Warburton (ref. 15) using an exact solution. 

(2) An orthotropic cylinder with freely supported edges investigated by Hoppmann 
(ref. 16) using an exact solution. 

(3) An isotropic 120° conical frustum with both free-free and clamped-free edges, 
investigated by Naumann (ref. 1) using a Rayleigh-Ritz analysis. 

(4) An isotropic shell having positive Gaussian curvature with freely supported edges 
investigated by Cooper (ref. 4) using a finite difference solution. 

(5) An isotropic shell having negative Gaussian curvature with freely supported edges 
also investigated by Cooper (ref. 4). 

(6) An isotropic annular plate with free edges, investigated by Raju (ref. 17) using 
an exact solution. 


Correlation With Previous Investigations 

Frequencies and mode shapes were computed for these shells, and the frequencies 
were compared with those from existing solutions. Ten elements were used to represent 
each shell. One hundred integration intervals within each element were used. In the cor- 
relations to follow, the quantities to be compared are called frequency parameters. For 
the cylinders, conical frustums, and annular plate, the parameter is the square of the 
circular frequency. For the shells of positive and negative Gaussian curvature, the param- 
eter is a dimensionless frequency defined in the appropriate tables and figures. 

Frequencies based on the methods of previous investigators were obtained as follows: 

(1) For the cylinders and the annular plate, the methods of the previous investigators 
(refs. 15, 16, and 17) were automated for computation on a digital computer. Some of the 
physical data of the orthotropic cylinder were obtained from reference 18. 

(2) For the conical frustums, the computer program of Naumann (ref. 1) was used. 

(3) For the shells of positive and negative Gaussian curvature, frequencies were 
provided by Paul A. Cooper who obtained them by use of a computer program based on the 
procedure described in reference 4. 
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Presentation of Results 


Calculations were performed with the present method for each shell using a range 
on the circumferential wave number n of 0 through 10 because the minimum frequencies 
for the shells were in this range. Since all calculations by the present method were based 
on representations by 10 elements, approximately 60 modes per value of n were gener- 
ated. Because presentation of all these modes appeared impractical, a decision had to be 
made as to which modes to present for each shell. A main consideration in selection of 
a mode for presentation was whether a calculation of the frequency parameter was avail- 
able from one of the methods of the previous investigators for comparison. 

Minimum frequencies .- The minimum frequency parameter for each value of n 
was available for every shell. Correlations for minimum frequency parameters are pre- 
sented in table IV. The information in this table is shown graphically in figures 5 to 11. 
The frequency parameters in these figures should be viewed both as results of the present 
analysis and as results of the previous investigators since the differences are too small 
to be seen on the plots. Some experimental results from references 16 and 1 are shown 
in figures 6 and 7, respectively. 

Higher modes of cylinders .- The selection of which higher frequency parameters to 
present was made for each shell on an individual basis. For the cylinders, the solutions 
of Arnold and Warburton and of Hoppmann give exact frequency parameters for all modes. 
In order to describe which of these modes were selected for correlation, it is necessary 
to discuss briefly the nature of the exact mode shapes for the freely supported cylinders. 
The exact mode shapes have the form: 

u = A mn cos cos n0 
L 

v = B mn sin S i n n 0 

J_i 

w = C mn sin cos nd 
L 

where m takes on integer values and A mn , B mn , and C mn are constants which char- 
acterize a mode. 

Thus, m is equal to the number of nodal circles in u or one plus the number of 
nodal circles in v or w. For a given pair of values for m and n, three modes are 
possible. Each mode corresponds to a different ratio of A mn : B m n • Cmn- For each 
value of n, modes were arbitrarily selected for values of m from 1 to 5. The cor- 
responding frequency parameters computed by the present analysis were identified by 
inspecting the computed mode shapes and counting nodal circles. Correlations for the 
higher frequency parameters of the cylinders are presented in tables V and VI for the 
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values of n and m considered. For the purpose of illustration, the variations with m 
of the three frequency parameters associated with each value of m are shown in fig- 
ures 12 and 13 for n = 2. 

Higher modes of conical frustums .- For the conical frustums, a large number of 
frequency parameters were available from Naumann's Rayleigh- Ritz procedure for com- 
parison with results of the present analysis. A correlation was made for n = 2. Modes 
were selected for presentation as follows: First, all modes having from one to five nodal 

circles in the w-displacement were examined. For the freely supported cylinders, there 
were exactly three modes having a given number of circumferential waves and a given 
number of circular node lines in the w-displacement. For the conical frustums, there 
are usually not exactly three. For purposes of presenting frequencies in the present 
analysis, an arbitrary selection of modes was made. If there are three or less modes 
corresponding to a given number of nodal circles in w, all are presented. If there are 
more than three, only three are presented. To avoid any misunderstanding, the following 
table tells which modes of the conical frustums having five or less nodal circles in w 
are not presented: 



Nodal circles 
in w 

Nodal circles 
in u 

Nodal circles 
in v 

Free-free 

1 

0 

0 

Clamped-free 

0 

1 

0 


3 

0 

1 


4 

2 

5 


3 

4 

7 


4 

6 

8 


The correlations for the higher frequency parameters of the conical frustums for n = 2 
are shown in tables VII and VIII. This information is also shown graphically in figures 14 
and 15. As in figures 5 to 11, the results in figures 14 and 15 can be interpreted as either 
the results of the present analysis or the other investigation since the results are coinci- 
dent for plotting purposes. As a matter of interest, the mode shapes from the present 
analysis corresponding to the frequency parameters in figures 14 and 15 are shown in fig- 
ures 16 and 17. 

Higher modes of shells having positive and negative Gaussian curvature .- For the 
shells of positive and negative Gaussian curvature, no higher frequency parameters were 
available from Cooper's method for correlation. As a matter of interest, some higher 
frequency parameters computed by the present analysis for these shells are presented for 
n = 2 in tables IX and X. Specifically, frequency parameters for modes having four or 
less circular node lines in w are presented. As was the case with the cylinders, there 
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are exactly three frequencies corresponding to each number of circular node lines in w. 
The variation of the three frequency parameters with the number of circular node lines 
in w is shown in figures 18 and 19. The mode shapes corresponding to the frequency 
parameters in figures 18 and 19 are shown in figures 20 and 21. 

Discussion of the Applications 

The correlations in tables IV to VIII can be summarized as follows: All frequency 
parameters presented for the two cylinders, the shells of positive and negative Gaussian 
curvature, and the annular plate showed agreement with the results of the previous inves- 
tigators through at least the second significant figure. For the cylinders, the agreement 
was in most cases through six significant figures. For the free-free conical frustum, 
most frequency parameters predicted by the present analysis agreed with the frequency 
parameters from the method of reference 1 through at least the second significant figure. 
The exceptions were the ninth, eleventh, and twelfth values of u> 2 listed in table VII in 
which the frequency parameters from the present analysis were lower in the second sig- 
nificant figure. It is noted in reference 19 that under certain conditions (that are met by 
the present analysis) the finite- element method is equivalent to the Rayleigh- Ritz method 
in that both methods give upper bounds to the exact frequencies. Therefore, it follows 
that the frequency parameters predicted by the present analysis are better approxima- 
tions to the corresponding exact frequency parameters than are the frequency parameters 
from the method of reference 1. It is believed that the first significant figure in the three 
frequency parameters is probably correct and that the lack of agreement for these fre- 
quency parameters does not indicate any significant inaccuracy in the results of the pres- 
ent analysis. For the clamped-free conical frustum, most frequency parameters pre- 
dicted by the present analysis again agreed with the frequency parameters from the 
method of reference 1 through at least the second significant figure. The exceptions in 
this instance are the minimum frequency parameters in table IV for n = 1, 2, and 3 and 
the eighth and ninth values of the frequency parameters o> 2 listed in table VIII. For the 
eighth and ninth values of u> 2 listed in table VIII, the present analysis predicted fre- 
quency parameters which were lower than the corresponding values from the method of 
reference 1. As with the free-free conical frustum, these differences occurred in the 
second significant figure. It is again concluded that the present results are closer to the 
exact frequency parameters. In the case of the noted disagreement in table IV, the fre- 
quency parameters from the present analysis are higher (in the first significant figure 
for n = 2) than the corresponding result from the method of reference 1. By reasoning 
similar to that used in the preceding discussion, it is concluded that the present analysis 
is somewhat inaccurate for the minimum frequency parameters of the clamped-free coni- 
cal frustum for n = 1, 2, and 3. 
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No dependence of accuracy on the circumferential wave number was noted. Further- 
more, only slight degradation of accuracy was noted as modes with increasing numbers of 
circular node lines in w were considered. In the latter instance, significant degrada- 
tion might have been expected. The only notable degradation of accuracy was with regard 
to edge constraint. The frequency parameters of the conical frustum with a clamped edge 
(tables IV and VIII) showed worse correlation with the results of the previous investiga- 
tions than did the frequency parameters of shells with free or freely supported edges. 

For the reasons stated in the introduction, no correlations are presented for mode 
shapes. However, a cursory correlation between the mode shapes from the present anal- 
ysis and such mode shapes as were available from the methods of the previous investiga- 
tors was made. The computed mode shapes for the cylinders appeared to coincide with 
the exact mode shapes which are sine and cosine curves. The computed mode shapes 
corresponding to some of the minimum frequencies of the free-free conical frustum 
appeared to agree with those mode shapes published in reference 1. The authors also 
made some correlations of mode shapes for a few of the higher modes of the conical 
frustums obtained by the method of reference 1. For these mode shapes the present anal- 
ysis and the method of reference 1 appeared to agree very well. 

Computational Efficiency and Reliability 

As was stated in the introduction, two major objectives of the computer program 
were: (1) machine efficiency, that is, the ability to compute quickly a large number of 

frequencies and mode shapes; (2) reliability, that is, capability of predicting every mode 
in the range of the frequency spectrum of interest. Machine efficiency was achieved. 
Typically, over 600 frequencies and modal columns are computed in less than 15 minutes 
on the Control Data 6600 computer system. Reliability was not proven but is indicated by 
the correlations obtained by the exact theory for the cylinders. The eigenvalue problems 
generated were well conditioned, since only single-precision arithmetic was required for 
accurate solution. 


Limitations 

Experience with this computer program and with the correlations are believed to 
indicate that the major limitation of the program in its present form lies in the approxi- 
mation of the normal displacement w by a third-order curve over each element. (See 
eq. (7).) From the assumption that w is a third-order curve in each element, it follows 
that: 

(1) Certain moment resultants are discontinuous across junctures between elements. 
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(2) The moment distributions may be seriously in error if the moment distribution 
in a region represented by a single element has a variation of higher order than linear. 

Some results of modal stress calculations for modes with large moment gradients 
are suspected to be in error. Furthermore, the errors previously noted in some mini- 
mum frequencies of the clamped-free conical frustum are believed to stem from the third- 
order approximation since a steep moment gradient is known to occur near a clamped 
edge. A possible remedy is to increase by two the order of polynomial representation of 
w and to require continuity of curvature across element junctures. 

Another limitation which is emphasized is the restriction of the present analysis to 
shells for which the shell surface does not intersect the axis of the shell. Thus, this 
analysis is not applicable to configurations such as a hemisphere. 

Finally, the reader is reminded that the analysis is restricted to shells with con- 
tinuous stiffness distributions as noted from the conditions imposed by equation (33). 

This restriction is easily removed by replacing the last two equalities of equation (33) by 
appropriate conditions on the continuity of stress and moment resultants across element 
junctures. 


CONCLUDING REMARKS 

An analytical procedure based on the finite-element method is developed for com- 
puting natural frequencies and mode shapes of thin shells of revolution. The shells may 
have general meridional curvature and orthotropic elastic properties. The details of a 
computer program based on this procedure are described. 

A distinguishing feature of the procedure is that it employs an element which is 
geometrically exact in that the actual geometry of the shell being analyzed is input to the 
analysis in the form of functions. The displacements of the shell within an element are 
approximated by third-order polynomials which are defined over the element. Inter- 
element compatibility is expressed by equating displacements and rotations at all junc- 
tures between elements. The required integrations for computing the element stiffness 
and mass matrices are performed numerically by using the trapezoidal rule. The stiff- 
ness and mass matrices for the complete shell are formed by superposition. Edge con- 
straints are incorporated by deleting rows and columns from the complete shell stiffness 
and mass matrices. The resulting symmetric eigenvalue problem is solved by a standard 
method. 
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The computer program has been applied to several shells: 

(1) An isotropic cylinder with freely supported edges 

(2) An orthotropic cylinder with freely supported edges 


(3) A 120° conical frustum with free-free edges 

(4) A 120° conical frustum with clamped-free edges 

(5) A shell having positive Gaussian curvature with freely supported edges 

(6) A shell having negative Gaussian curvature with freely supported edges 

(7) An annular plate with free-free edges 

The main results and conclusions are as follows: 

1. Very generally, excellent agreement was noted between frequencies from the 
present analysis and frequencies from the previous investigations. 

2. The only inaccuracies of the present analysis which might be considered signifi- 
cant occurred in three minimum frequencies of the clamped-free conical frustum. This 
inaccuracy is believed to stem from the inability of third-order polynomials to conform 
to a steep stress gradient near a clamped edge, and consequently increasing the repre- 
sentation of the normal-displacement component to a fifth-order polynomial would be 
expected to result in overall excellent agreement. 

3. The computer program performs with very short running times and no modes are 
overlooked in computation. 

4. The natural frequencies and mode shapes from this method appear to constitute 
reliable input for forced response calculations for structures involving shells of 
revolution. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 3, 1968, 

124-08-05-08-23. 
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APPENDIX 


ELEMENTS OF MATRIX [r] 

(See eq. (16)] 

The elements of matrix [r] are as follows: 

_ C n r C 12 r C 22 r R 22 n4 D 66 n2 ( r ’) 2 2K 12 n2 2K 22 n2 

11 R x 2 + r 1 r 2 + R 2 2 r 3 r 3 + rR x + rR 2 


(Al) 


r 12 = r 21 = " 


° 22 r ’ n2 D 66 n2r ' K 12 r ’ K 22 r ' 


R 1 r 2 


(A2) 


^ D °12 n2 K ll r K 12 r 

r 13 = r 31 r 


r 14 = r 41 = 


R 1 r 2 

Ci 2 r' C 22 r’ D 12 n 2 R] D 22 r’n 2 DggnV K n R\r K 12 r’ 


R i Rn 


rRi 2 r 2 Ri r 2 Ri 


Ri 


R 2 
K 1 


(A3) 


K 12 rR l , K 22 n2r ’ , K 22 r ’ , K 66 n2r ’ 
R x 2r 2 r 2 r 1 r 2 + r 2 


(A4) 


D D C ll r . C 12 r . D 12 n2 . K ll r . K 12 n2 . K 12 r 

R i * = R *i = — — + — — + — - — + 9- + — f — + r~ r" 

R x 2 r 1 r 2 


1.1 c = Kci = — + — + 

15 o-L Rj R 2 rR^ 


(A 5) 


D _ n _ Ci 2 n C 22 n D 22 n 3 Dggn(r') 2 K 12 n K 22 n 3 K 22 n 

R 16 = R 6 1 + -Rf + - - 72 ^- + RlRi + ~ + i^ 2 - 

K 6 6"( r '> 2 


(A6) 


R 22 " 


R71 rR 2 r 

(A7) 

D 22( r ’) 2 °66 n2 
r 1 r 

(A8) 

R 32 = D 12 r ’ 

(A9) 
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D 12 R l r ' D 22( r ’) 2 ° 66 n * K 22 (r')^ K 66 n^ 

R 24 42 „ 2 " rRi ~ rRi r ’ r 


,t\2 


Ri 


r 25 “ R 52 - ” 


Dl 2 r ’ 

R 1 


- K 12 r * 


D 22 nr ’ D 66 rn K 22 nr ’ K 66 nr ’ 


r 26 - R 62 - " rR2 rR2 r r 


D 66 n 


r 27 = R 72 ~ + K 66 n 


R 33 “ D ll r 


Ri i R i r 

Rq/j = R/i q = — “ Ki o r ’ " 


D 12 r ' 


34 " “43 - " 12 Rj 


Di 1 T 

R 35 = R 53 = - — - K - r 


R, 


k ll 


R _ R _ ° 12 n _ K n 
R 36 " R 63 " “r^" K 12 n 


R 37 _ R 73 " 0 


R. 


C 22 (-") 2 , C 6 6° 2 , D ll( R i) 2r 2D 12 R i r ' , P 22( r ’) 2 , P 66" 


44 r 


R l 4 


R ^ 2 rR ^ 2 rRj 2 


2 K 12 r ' R l , 2 K 22< r '> 2 , 2 K 66 n2 
R x 2 + rR l rR l 


r 45 " R 54 _ C 12 r 


, _ D n R i r { D l2 r ' K ll R l r [ 2K 12 r ' 


Rt 6 Ri 


Ri 


2 Ri 


(A10) 

(All) 

(A 12) 

(A 13) 
(A14) 

(A 15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 
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- r _ C 22 r n , C 66 nr ’ D 12 nR l D 22 nr ’ D 66 r ’ n K 12 nR l 
46 ' 64 " r r " R^Rg + rR i R 2 + rR l R 2 " R x 2 

| K 22 nr ’ [ K 22 nr ’ [ K 66 nr ' t K 66 nr ' 
rRj rR 2 + rRj + r R 2 


R 47 “ R 74 ” “ C 66 n " 


p 66 n K 66 n _ K 66 n 

RjR2 Rj R 2 


R -C r + D H r + 2K ll r 
R 55 “ C n r + TTT + 


R, 


Rt 


R - R -Cm ° 12n . Kl2n . Kl2n 
R 56 - R 65 - C 12 n + + + 


R 57 _ R 75 ~ 0 


(A21) 

(A22) 

(A23) 

(A24) 

(A25) 


_ C 22 n2 C 66( r ’) 2 D 22 n2 D 66( r ’) 2 2 K 22° 2 2 I W r ') 2 


R 66 r + r 


rR 2 2 rR 2 2 rR 2 


rR<] 


D fiR r’ 2K RR r’ 
r 67 = r 76 - -Cee 1 " - ~2 — 


Rc] 


Rc 


R - C r 4. ° 66r 4. 2K 66 r 
r 77 " c 66 r + tty 


Rri 


r 2 


(A26) 


(A27) 


(A28) 
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TABLE III.- EDGE CONSTRAINTS 


ICASE 

Description 

Equations for edge constraint 

Rows and columns deleted 

1 

Free-free 

None 

None 

2 

Free — freely supported 

v(L) = w(L) = 0 

(6K + 1), (6K + 3) 

3 

Freely supported — free 

v(0) = w(0) = 0 

1, 3 

4 

Free— simply supported 

u(L) = v(L) = w(L) = 0 

(6K + 1), (6K + 2), (6K + 3) 

5 

Simply supported— free 

u(0) = v(0) = cu(0) = 0 

1, 2, 3 

6 

Free- clamped 

u(L) = v(L) = w(L) = /3(L) = 0 

(6K + 1), (6K + 2), (6K + 3), (6K + 4) 

7 

Clamped-free 

u(0) = v(0) = w(0) = /3(0) = 0 

1, 2, 3, 4 

8 

Freely supported— 
freely supported 

v(0) = w(0) = 0 
v(L) = w(L) = 0 

1, 3, (6K + 1), (6K + 3) 

9 

Simply supported— 
simply supported 

u(0) = v(0) = w(0) = 0 
u(L) = v(L) = w(L) = 0 

1, 2, 3, (6K + 1), (6K + 2), (6K + 3) 

10 

Clamped-clamped 

u(0) = v(0) = w(0) = /3(0) = 0 
u(L) = v(L) = w(L) = /3(L) = 0 

1, 2, 3, 4, (6K + 1), (6K + 2), (6K + 3), (6K + 4) 

11 

Freely supported— 
simply supported 

v(0) = w(0) = 0 
u(L) = v(L) = w(L) = 0 

1, 3, (6K + 1), (6K + 2), (6K + 3) 

12 

Freely supported— 
clamped 

v(0) = w(0) = 0 
u(L) = v(L) = w(L) = /3(L) = 0 

1, 3, (6K + 1), (6K + 2), (6K + 3), (6K + 4) 

13 

Simply supported— 
freely supported 

u(0) = v(0) = w(0) = 0 
v(L) = w(L) = 0 

1, 2, 3, (6K + 1), (6K + 3) 

14 

Simply supported— 
clamped 

u(0) = v(0) = w(0) = 0 
u(L) = v(L) = w(L) = j3(L) = 0 

1, 2, 3, (6K + 1), (6K + 2), (6K + 3), (6K + 4) 

15 

Clamped — 

freely supported 

u(0) = v(0) = w(0) = /3(0) = 0 
v(L) = w(L) = 0 

1, 2, 3, 4, (6K + 1), (6K + 3) 

16 

Clamped— 

simply supported 

u(0) = v(0) = w(0) = /3(0) = 0 
u(L) = v(L) = w(L) = 0 

1, 2, 3, 4, (6K + 1), (6K + 2), (6K + 3) 
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TABLE IV.- COMPARISON OF MINIMUM FREQUENCY PARAMETERS AS COMPUTED BY 
PRESENT METHODS WITH THOSE OF PREVIOUS INVESTIGATIONS 



sec~2 

Circumferential 
wave number, 
n 

Freely supported— freely supported cylinder 

Isotropic 

Orthotropic 


Present analysis 

Reference 15 

Present analysis 

Reference 16 

0 

3.86111 X 10 8 

3.86111 x 10 8 

8.01167 X 10 8 

8.01167 X 10 8 

1 

1.17339 X 10 8 

1.17339 X 10 8 

1.09559 X 10 8 

1.09559 X 10 8 

2 

2.25430 X 10 7 

2.25430 X 10 7 

2.15504 X 10 7 

2.15504 X 10 7 

3 

5.95827 X 10 7 

5.95827 X 10 7 

5.53142 X 10 7 

5.53142 X 10 7 

4 

2.17401 X lO 8 

2.17401 X lO 8 

1.85439 X 10 8 

1.85439 x 10 8 

5 

1.11765 X 10 6 

1.11765 X 10 6 

4.74829 X 10 8 

4.74829 X 10 8 

6 

9.09145 X 10 5 

9.09145 X 10 5 

1.01231 X 10 9 

1.01231 x 10 9 

7 

1.11505 X 10 6 

1.11505 X 10 6 

1.90756 X 10 9 

1.90756 X 10 9 

8 

1.64300 X 10 6 

1.64300 X 10 6 

3.29043 X 10 9 

3.29043 x 10 9 

9 

2.50514 X 10 6 

2.50514 X 10 6 

5.31084 X 10 9 

5.31084 X 10 9 

10 

3.75508 X 10 6 

3.75508 X 10 6 

8.13868 X 10 9 

8.13868 x 10 9 


o;2, sec-2 

Circumferential 
wave number, 
n 

120° conical frustum 

Free-free 

Clamped-free 


Present analysis 

Reference 1 

Present analysis 

Reference 1 

0 

0 

0 

2.5380 X 10 4 

2.5378 X 10 4 

1 

0 

0 

3.7049 x 10 5 

3.6125 X 10 5 

2 

2.8727 X 10 2 

2.8725 X 10 2 

7.2558 X 10 4 

6.9075 x 10 4 

3 

1.9154 X 10 3 

1.9149 X 10 3 

2.2666 X 10 4 

2.1638 X 10 4 

4 

6.3759 X 10 3 

6.3728 X 10 3 

1.3694 X 10 4 

1.3434 X 10 4 

5 

1.5047 X 10 4 

1.5038 X 10 4 

1.7636 X 10 4 

1.7611 X 10 4 

6 

2.7824 X 10 4 

2.7815 x 10 4 

2.8529 x 10 4 

2.8520 x 10 4 

7 

4.4394 X 10 4 

4.4387 X 10 4 

4.4524 x 10 4 

4.4507 X 10 4 

8 

6.6322 X 10 4 

6.6310 X 10 4 

6.6338 X 10 4 

6.6324 X 10 4 

9 

9.5417 X 10 4 

9.5394 X 10 4 

9.5418 X 10 4 

9.539 X 10 4 

10 

1.3333 X 10 5 

1.3329 X 10 5 

1.3333 X 10 5 

1.3329 X 10 5 
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TABLE IV.- COMPARISON OF MINIMUM FREQUENCY PARAMETERS AS COMPUTED BY 
PRESENT METHODS WITH THOSE OF PREVIOUS INVESTIGATIONS - Concluded 



Dimensionless frequency parameter, £2 = wR[(p/E)(l 

- fX 2 )] V2 

Circumferential 
wave number, 
n 

Shell of positive 
Gaussian curvature, 
freely supported edges 

Shell of negative 
Gaussian curvature, 
freely supported edges 


Present analysis 

Reference 4 

Present analysis 

Reference 4 

0 

— 

— 

0.640 

0.640 

1 

0.411 

0.412 

.368 

.368 

2 

.360 

.362 

.157 

.157 

3 

.340 

.340 

.0628 

.0628 

4 

.331 

.331 

.01970 

.01972 

5 

.327 

.327 

.00779 

.00784 

6 

.324 

.324 

.01923 

.01924 

7 

.323 

.322 

.02804 

.02805 

8 

.322 

.321 

.02580 

.02609 

9 

.321 

.321 

.0240 

.0241 

10 

.321 

.321 

.0292 

.0292 

Circumferential 
wave number, 

sec“2 

Plates 

n 

Present analysis 

Reference 17 

0 

86.74 

86.74 

1 

295.8 

295.8 

2 

18.24 

18.24 

3 

130. 

,5 

130.5 

4 

443.8 

443.8 

5 

1 087 


1 087 


6 

2 215 


2 215 


7 

4 003 


4 003 


8 

6 660 


6 660 


9 

10 415 


10 415 


10 

15 532 


15 532 
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TABLE V.- CORRELATION OF FREQUENCY PARAMETER w 2 OF FREELY SUPPORTED CYLINDER 
AS COMPUTED BY PRESENT ANALYSIS AND EXACT SOLUTION OF REFERENCE 15 


a) 2 , sec" 


Reference 15 

Present analysis 

Reference 15 

Present analysis 

Reference 15 

0 



n = 

1 





n = : 

l 


1.950193 x 

10 9 

1.17339 

X 

10 8 

1.17339 

X 

10 8 

2.25430 

X 

10 7 

2.25430 x 

10 7 

9.00656 

X 

10 8 

1.41017 

X 

10 9 

1.41017 

X 

10 9 

3.1160 

X 

10 9 

3.116 x 

10 9 

3.86111 

X 

10 8 

4.09860 

X 

10 9 

4.09860 

X 

10 9 

9.65443 

X 

10 9 

9.65443 x 

10 9 

1.53415 

X 

10 9 

5.83357 

X 

10 8 

5.83356 

X 

10 8 

1.90073 

X 

10 8 

1.90073 x 

10 8 

4.47961 

X 

10 9 

2.78500 

X 

10 9 

2.78500 

X 

10 9 

4.59200 

X 

10 9 

4.59200 X 

10 9 

1.54445 

X 

10 9 

6.67902 

X 

10 9 

6.67902 

X 

10 9 

1.24328 

X 

10“ 

1.24328 x 

10™ 

1.58626 

X 

10 9 

1.02452 

X 

10 9 

1.02451 

X 

10 9 

4.75299 

X 

10 8 

4.75290 x 

10 8 

9.96570 

X 

10 9 

1.18772 

X 

10“ 

1.18772 

X 

10“ 

6.60278 

X 

10 9 

6.60268 X 

10 9 

3.47500 

X 

10 9 

4.51464 

X 

10 9 

4.51457 

X 

10 9 

1.75057 

X 

10“ 

1.75056 X 

10 10 

1.59995 

X 

10 9 

1.26785 

X 

10 9 

1.26782 

X 

10 9 

7.55161 

X 

10 8 

7.55117 x 

10 8 

1.75655 

X 

10 8 

1.94100 

X 

10“ 

1.94088 

X 

10 10 

2.49175 

X 

10“ 

2.49164 X 

10 10 

6.17778 

X 

10 9 

7.05640 

X 

10 9 

7.95582 

X 

10 9 

9.22905 

X 

10 9 

9.22834 X 

10 9 

1.60571 

X 

10 9 

1.39208 

X 

10 9 

1.39201 

X 

10 9 

9.73034 

X 

10 8 

9.72909 x 

10 8 

2.73485 

X 

10 10 

2.91681 

X 

10“ 

2.91617 

X 

10“ 

3.46020 

X 

10“ 

3.45961 X 

10™ 

9.65279 

X 

10 9 

1.04453 

X 

10 “ 

1.04425 

X 

10“ 

1.25978 

X 

10“ 

1.25946 X 

10™ 

3 



n = 

4 



n = 5 

5.95827 

X 

10 6 

2.17401 

X 

10 6 

2.17401 

X 

10 6 

1.11765 

X 

10 6 

1. 11765 X 

10 6 

6.15055 

X 

10 9 

1.04833 

X 

10“ 

1.04833 

X 

io“ 

1.60902 

X 

10 10 

1.60902 X 

io“ 

1.85829 

X 

10 10 

3.09779 

X 

10“ 

3.09779 

X 

10“ 

4.68747 

X 

10™ 

4.68747 x 

10“ 

6.66563 

X 

10 7 

2.70603 

X 

10 7 

2.70601 

X 

10 7 

1.27354 

X 

10 7 

1.27352 x 

10 7 

7.56576 

X 

10 9 

1.18310 

X 

10“ 

1.18310 

X 

10“ 

1.73881 

X 

10 10 

1.73881 x 

10“ 

2.15282 

X 

1010 

3.40266 

X 

10“ 

3.40266 

X 

10“ 

4.99865 

X 

1010 

4.99865 X 

10“ 

2.13118 

X 

10 8 

1.00599 

X 

10 8 

1.00594 

X 

10 8 

5.13733 

X 

10 7 

5.13704 X 

10 7 

9.67963 

X 

10 9 

1.39500 

X 

10 10 

1.39499 

X 

10“ 

1.94829 

X 

10 10 

1.94828 x 

10“ 

2.66367 

X 

10“ 

3.92032 

X 

10“ 

3.92031 

X 

10“ 

5.52220 

X 

10*0 

5.52219 x 

10“ 

4.10708 

X 

10 8 

2.22493 

X 

10 8 

2.22462 

X 

10 8 

1.24623 

X 

10 8 

1.24601 x 

10 8 

1.24372 

X 

10 10 

1.67738 

X 

10 10 

1.67733 

X 

10“ 

2.23239 

X 

10 10 

2.23230 x 

10“ 

3.39978 

X 

10 10 

4.65653 

X 

10“ 

4.65745 

X 

io“ 

6.26256 

X 

10™ 

6.26249 x 

10 10 

6.12119 

X 

10 8 

3.72329 

X 

10 8 

3.72208 

X 

10 8 

2.27379 

X 

10 8 

2.27282 X 

10 8 

4.36131 

X 

10 10 

2.02963 

X 

10“ 

2.02922 

X 

10“ 

7.22276 

X 

10 10 

7.22237 x 

10“ 

1.58813 

X 

10“ 

5.61740 

X 

10“ 

5.61693 

X 

10“ 

2.58897 

X 

10™ 

2.58853 x 

io“ 


Present analysis 


1.950193 x 10 9 
10 8 
10 8 

10 9 
10 9 
10 9 

10 9 
10 9 


9.00656 

3.86111 

1.53415 

4.57961 

1.54445 

1.58626 

9.96586 

3.47507 

1.59995 

1.75668 

6.17832 

1.60572 

2.73552 

9.65539 


10 


10 " 

10 8 

10 9 

in9 


10 


5.95827 

6.15055 

1.85829 

6.66567 

7.56577 

2.15282 

2.13124 

9.67973 

2.66367 

4.10749 

1.24380 

3.39988 

6.12256 

4.36214 

1.58850 


x 10 6 
X 10 9 
x 10 10 

X 10 7 
X 10 9 
x 10 10 

X 10 8 
x 10 9 
x 10 10 

X 10 8 


x 10 


10 


x 10 10 

x 10 8 
x 10 10 
x 10 10 
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TABLE V.- CORRELATION OF FREQUENCY PARAMETER to 2 OF FREELY SUPPORTED CYLINDER 
AS COMPUTED BY PRESENT ANALYSIS AND EXACT SOLUTION OF REFERENCE 15 - Concluded 



a) 2 , sec -2 

m 

Present analysis 

Reference 15 

Present analysis 

Reference 15 

Present analysis 

Reference 15 


n = 

6 

n = 

7 

n = 

8 


r 9.09145 X 10 5 

9.09145 x 10 5 

1.11505 x 10 6 

1.11505 X 10 6 

1.64300 x 10 6 

1.64300 x 10 6 

1 

< 2.29598 x 10 10 

2.29598 x 10 10 

3.10870 X 10 10 

3.10870 X 10 10 

4.04694 X 10 10 

4.04694 x 10 10 


V. 6.62863 X 10 10 

6.62863 x 10 10 

8.92182 x 10 10 

8.92182 X 10 10 

1.15673 x 10 11 

1.15673 x 10 11 


r 6.96234 x 10 6 

6.96226 x 10 6 

4.56899 X 10 6 

4.56892 X 10 6 

3.77396 X 10 6 

3.77391 x 10 6 

2 

< 2.42233 x 10 10 

2.42233 X 10 10 

3.23267 X 10 10 

3.23267 X 10 10 

4.16922 X 10 10 

4.16922 X 10 10 


L 6.94380 x 10 10 

6.94380 x 10 10 

9.23964 X 10 10 

9.23964 X 10 10 

1.18870 x 10 11 

1.18870 X 10 11 


C 2.85635 X 10 7 

2.85615 X 10 7 

1.74106 X 10 7 

1.74092 X 10 7 

1.18524 X 10 7 

1.18514 x 10 7 

3 

< 2.62905 X 10 10 

2.62903 X 10 10 

3.43697 X 10 10 

3.43696 X 10 10 

4.37157 X 10 10 

4.37156 X 10 10 


7.47182 x 10 10 

7.47182 x 10 10 

9.77096 X 10 10 

9.77096 X 10 10 

1.24207 x 10 11 

1.24207 x 10 11 


C 7.33039 x 10 7 

7.32876 X 10 7 

4.56222 x 10 7 

4.56103 X 10 7 

3.03024 x 10 7 

3.02935 x 10 7 

4 

< 2.91264 x 10 10 

2.91254 x 10 10 

3.71926 X 10 10 

3.71916 x 10 10 

4.65240 X 10 10 

4.65231 X 10 10 


L 8.2 1 554 x 10 10 

8.21548 X 10 10 

1.05176 x 10 11 

1.05176 x 10 11 

1.31698 x 10 11 

1.31698 x 10 11 


C 1.42224 x 10 8 

1.42150 x 10 8 

9.20757 x 10 7 

9.20184 x 10 7 

6.21633 X 10 7 

6.21190 x 10 7 

5 

< 3.27094 x 10 10 

3.27049 x 10 10 

4.07780 X 10 10 

4.07733 x 10 10 

5.01043 X 10 10 

5.00996 X 10 10 


L 9.17738 x 10 10 

9.17704 x 10 10 

1.14813 X 10 11 

1.14812 X 10 11 

1.41378 X 10 11 

1.41353 X 10 11 

m 

n = 

9 

n = 

10 



r 2.50514 x 10 6 

2.50514 x 10 6 

3.75508 X 10 6 

3.75508 X 10 6 



1 

< 5.11057 X 10 10 

5.11057 X 10 10 

6.29952 x 10 10 

6.29952 x 10 10 




L 1.45653 X 10 11 

1.45653 x 10 11 

1.79157 X 10 11 

1.79157 x 10 11 




f 3.91728 x 10 6 

3.91725 x 10 6 

4.75813 X 10 6 

4.75810 x 10 6 



2 

< 5.23161 X 10 10 

5.23162 X 10 10 

6.41965 x 10 10 

6.41965 X 10 10 




L 1.48862 X 10 11 

1.48862 x 10 11 

1.82376 x 10 11 

1.82376 X 10 11 




C 9.27100 X 10 6 

9.27021 x 10 6 

8.48 546 X 10 6 

8.48484 X 10 6 



3 

< 5.43242 x 10 10 

5.43241 X 10 10 

6.61922 X 10 10 

6.61921 X 10 10 




L 1.54218 X 10 11 

1.54218 X 10 11 

1.87746 x 10 11 

1.87746 x 10 11 




f 2.17922 X 10 7 

2.17854 x 10 7 

1.72993 x 10 7 

1.72939 X 10 7 



4 

< 5.71189 x 10 10 

5.71180 X 10 10 

6.89750 X 10 10 

6.89741 X 10 10 




V. 1.61727 x 10 11 

1.61727 x 10 11 

1.95275 X 10 11 

1.95272 x 10 11 




r 4.41539 x 10 7 

4.41191 X 10 7 

3.34003 X 10 7 

3.33725 x 10 7 



5 

< 6.06912 X 10 10 

6.06864 x 10 10 

7.25384 X 10 10 

7.25336 x 10 10 




L 1.71401 X 10 11 

1.71401 x 10 11 

2.04966 x 10 11 

2.04960 x 10 11 
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TABLE VI.- COMPARISON OF FREQUENCY PARAMETER u 2 OF AN ORTHOTROPIC CYLINDER 
WITH FREELY SUPPORTED EDGES AS COMPUTED BY PRESENT ANALYSIS 
AND EXACT METHOD OF REFERENCE 16 



w 2 , sec~2 

m 

Present analysis 

Reference 16 

Present analysis 

Reference 16 

Present analysis 

Reference 16 


n = 

0 

n = 

: 1 

n = 

: 2 


C 8.37907 x 10 9 

8.37907 x 10 9 

1.09559 x 10 8 

1.09559 x 10 8 

2.15504 x 10 7 

2.15504 x 10 7 

1 

( 2.00989 x 10 9 

2.00989 x 10 9 

6.98986 x 10 9 

6.98986 x 10 9 

2.18774 x 10 10 

2.18774 x 10 10 


L 8.01167 x 10 8 

8.01167 x 10 8 

1.76461 x 10 10 

1.76461 x 10 10 

4.35227 X 10 10 

4.35227 x 10 10 


r 6.73909 x 10 9 

6.73909 X 10 9 

9.33628 x 10 8 

9.33624 x 10 8 

1.96173 x 10 8 

1.96170 x 10 8 

2 

< 1.00045 x 10 10 

1.00045 x 10 10 

1.05341 x 10 10 

1.05341 x 10 10 

2.39218 x 10 10 

2.39218 x 10 10 


L 3.20467 x 10 9 

3.20467 x 10 9 

2.20377 x 10 10 

2.20377 x 10 10 

5.00686 x 10 10 

5.00686 x 10 10 


r 7.77840 x 10 9 

7.77835 X 10 9 

2.47618 x 10 9 

2.47610 x 10 9 

7.07832 X 10 8 

7.07771 x 10 8 

3 

< 1.95747 X 10 10 

1.95744 x 10 10 

3.14843 x 10 10 

3.14841 x 10 10 

2.70286 x 10 10 

2.70284 x 10 10 


V, 7.21064 x 10 9 

7.21050 x 10 9 

1.41632 x 10 10 

1.41630 x 10 10 

6.10771 x 10 10 

6.10769 x 10 10 


r 7.97769 x 10 9 

7.97727 X 10 9 

4.14783 x 10 9 

4.14727 X 10 9 

1.56673 x 10 9 

1.56621 x 10 9 

4 

< 3.42687 x 10 10 

3.42661 x 10 10 

4.59040 x 10 10 

4.59018 x 10 10 

7.64494 x 10 10 

7.64474 X 10 10 


L 1.28198 x 10 10 

1.28198 x 10 10 

1.85782 x 10 10 

1.85770 x 10 10 

3.13158 x 10 10 

3.13143 x 10 10 


C 8.18535 x 10 9 

8.18303 x 10 9 

5.52799 x 10 9 

5.52532 X 10 9 

2.66133 X 10 9 

2.65858 x 10 9 

5 

< 5.32836 x 10 10 

5.32704 x 10 10 

6.48262 x 10 10 

6.48142 x 10 10 

9.60704 X 10 10 

9.60598 x 10 10 


L 2.00346 x 10 10 

2.00292 x 1010 

2.47182 X 10 10 

2.47122 X 10 10 

3.70574 X 10 10 

3.70507 x 10 10 

m 

n = 

3 

n = 

: 4 

n = 

5 


C 5.53142 x 10 7 

5.53142 x 10 7 

1.85439 x 10 8 

1.85439 x 10 8 

4.74829 x 10 8 

4.74829 x 10 8 

1 

< 4.77293 x 10 10 

4.77293 x 10 10 

8.43972 x 10 10 

8.43972 x 10 10 

1.31754 x 10 11 

1.31754 x 10 11 


L 8.54639 x 10 10 

8.54 639 x 10 10 

1.43693 X 10 11 

1.43693 X 10 11 

2.18346 X 10 11 

2.18346 x 10 11 


r 1.16635 X 10 8 

1.16633 x 10 8 

2.30036 X 10 8 

2.30034 x 10 8 

5.25903 X 10 8 

5.25901 x 10 8 

2 

< 4.85733 x 10 10 

4.85733 x 10 10 

8.44403 x 10 10 

8.44403 x 10 10 

1.31272 x 10 11 

1.31272 x 10 11 


L 9.33321 x 10 10 

9.33321 x 10 10 

1.52391 x 10 11 

1.52391 x 10 11 

2.27577 x 10 11 

2.27577 x 10 11 


C 3.18483 X 10 8 

3.18432 x 10 8 

3.55839 x 10 8 

3.55792 X 10 8 

6.44414 x 10 8 

6.44369 x 10 8 

3 

< 5.04809 X 10 10 

5.04807 x 10 10 

8.52479 x 10 10 

8.52477 x 10 10 

1.31219 x 10 11 

1.31219 x 10 11 


L 1.05863 x 10 11 

1.05864 x 10 11 

1.66119 x 10 11 

1.66119 x 10 11 

2.42199 x 10 11 

2.42198 x 10 11 


C 7.36447 x 10 8 

7.35971 x 10 8 

6.21337 x 10 8 

6.20889 x 10 8 

8.74824 x 10 8 

8.74385 x 10 8 

4 

< 5.37677 x 10 10 

5.37662 x 10 10 

8.73284 x 10 10 

8.73268 x 10 10 

1.32196 x 10 11 

1.32194 x 10 11 


L 1.22697 x 10 11 

1.22696 x 10 11 

1.84339 x 10 11 

1.84338 x 10 11 

2.61593 x 10 11 

2.61592 x 10 11 


C 1.40603 x 10 9 

1.40340 x 10 9 

1.08415 x 10 9 

1.08160 x 10 9 

1.27123 x 10 9 

1.26872 x 10 9 

5 

< 1.43650 x 10 11 

1.43641 x 10 11 

9.09817 x 10 10 

9.09739 x 10 10 

1.34611 x 10 11 

1.34603 x 10 11 


5.86279 x 10 10 

5.86202 x 10 10 

2.06739 X 10 11 

2.06733 X 10 11 

2.85407 x 10 11 

2.85339 x 10 11 
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TABLE VI.- COMPARISON OF FREQUENCY PARAMETER o> 2 OF AN ORTHOTROPIC CYLINDER 
WITH FREELY SUPPORTED EDGES AS COMPUTED BY PRESENT ANALYSIS 
AND EXACT METHOD OF REFERENCE 16 - Concluded 



u> 2 , sec"2 

m 

Present analysis 

Reference 16 

Present analysis 

Reference 16 

Present analysis 

Reference 16 


n = 

6 

n = 

7 

n = 

8 


C 1.01231 X 10 9 

1.01231 X 10 9 

1.90756 X 10 9 

1.90756 X 10 9 

3.29043 x 10 9 

3.29043 x 10 9 

1 

< 1.89737 X 10 11 

1.89737 X 10 11 

2.58318 x 10 11 

2.58318 x 10 11 

3.37482 x 10 11 

3.37482 x 10 11 


L 3.09484 X 10 11 

3.09484 X 10 11 

4.17137 x 10 11 

4.17137 X 10 11 

5.41321 X 10 11 

5.41321 x 10 11 


C 1.07881 X 10 9 

1.07880 X 10 9 

1.99482 x 10 9 

1.99481 x 10 9 

3.40265 x 10 9 

3.40265 x 10 9 

2 

< 1.88906 X 10 11 

1.88906 X 10 11 

2.47246 x 10 11 

2.47246 x 10 11 

3.36238 x 10 11 

3.36238 x 10 11 


L 3.19069 X 10 11 

3.19069 X 10 11 

4.26965 X 10 11 

4.26965 X 10 11 

5.51321 X 10 11 

5.51321 X 10 11 


r 1.21550 x 10 9 

1.215546 x 10 9 

2.16263 x 10 9 

2.16258 x 10 9 

3.61029 x 10 9 

3.61024 x 10 9 

3 

< 1.88205 X 10 11 

1.88205 X 10 11 

2.56061 X 10 11 

2.56061 X 10 11 

3.34688 x 10 11 

3.34688 x 10 11 


L 3.34351 X 10 11 

3.34351 X 10 11 

4.42737 x 10 11 

4.42737 X 10 11 

5.67465 x 10 11 

5.67462 X 10 11 


C 1.45929 X 10 9 

1.45885 X 10 9 

2.44390 x 10 9 

2.44346 X 10 9 

3.94411 x 10 9 

3.94365 X 10 9 

4 

< 1.88257 X 10 11 

1.88256 X 10 11 

2.55363 X 10 11 

2.55362 x 10 11 

3.33389 x 10 11 

3.33387 x 10 11 


L 3.54695 X 10 11 

3.54700 X 10 11 

4.63855 x 10 11 

4.63851 X 10 11 

5.89187 x 10 11 

5.89184 x 10 11 


C 1.85912 X 10 9 

1.85661 X 10 9 

2.88407 x 10 9 

2.88152 X 10 9 

4.44733 X 10 9 

4.44473 x 10 9 

5 

< 1.89540 X 10 11 

1.89533 x 10 11 

2.55664 X 10 11 

2.55657 x 10 11 

3.32856 X 10 11 

3.32849 x 10 11 


L 3.79646 X 10 11 

3.79631 x 10 11 

4.89809 x 10 11 

4.89791 x 10 11 

6.15984 x 10 11 

6.15972 x 10 11 

m 

n = 9 

n = 

10 



C 5.31084 X 10 9 

5.31084 X 10 9 

8.13868 x 10 9 

8.13868 X 10 9 



1 

< 4.27220 X 10 11 

4.27220 x 10 11 

5.27528 x 10 11 

5.27528 x 10 11 




L 6.82044 X 10 11 

6.82044 x 10 11 

8.39310 x 10 11 

8.39310 x 10 11 




C 5.45179 X 10 9 

5.45179 x 10 9 

8.31196 x 10 9 

8.31196 x 10 9 



2 

< 4.25852 X 10 11 

4.25852 X 10 11 

5.26066 X 10 11 

5.26066 X 10 11 




L 6.92169 X 10 11 

6.92169 x 10 11 

8.49529 x 10 11 

8.49529 x 10 11 




C 5.70636 X 10 9 

5.70631 x 10 9 

8.61988 x 10 9 

8.61982 X 10 9 



3 

( 4.24463 X 10 11 

4.24022 X 10 11 

5.24019 X 10 11 

5.24019 x 10 11 




L 7.08592 X 10 11 

7.08591 X 10 11 

8.66169 x 10 11 

8.66169 x 10 11 




C 6.10407 x 10 9 

6.10360 X 10 9 

9.09120 x 10 9 

9.09072 x 10 9 



4 

< 4.22240 X 10 11 

4.22239 x 10 11 

5.21849 X 10 11 

5.21849 x 10 11 




L 7.30801 X 10 11 

7.30800 X 10 11 

8.88772 x 10 11 

8.88769 X 10 11 




C 6.68682 X 10 9 

6.68415 X 10 9 

9.76702 x 10 9 

9.76427 x 10 9 



5 

l 4.21001 X 10 11 

4.21001 X 10 11 

5.20031 x 10 11 

5.20024 X 10 11 




L 7.58307 X 10 11 

7.58295 x 10 11 

9.16865 X 10 11 

9.16854 x 10 11 
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TABLE VII.- NATURAL FREQUENCIES OF A FREE- FREE 120° CONICAL 
FRUSTUM SHELL AS A FUNCTION OF NUMBER OF CIRCULAR 


NODE LINES IN NORMAL DISPLACEMENT w 

[„= 2 ] 


Number of circular 

Mode- identification 

sec 

-2 

node lines in w 

figure 

Present analysis 

Reference 1 

1 

16(a) 

3.903 x 10 4 

3.900 x 10 4 

1 

16(f) 

1.780 x 10 9 

1.781 x 10 9 

1 

16(j) 

4.337 x 10 9 

4.337 x 10 9 

2 

16(b) 

1.858 x 10 7 

1.861 X 10 7 

2 

16(k) 

8.045 X 10 9 

8.045 X 10 9 

3 

16(c) 

2.267 x 10 7 

2.286 x 10 7 

3 

16(g) 

2.422 x 10 9 

2.422 x 10 9 

3 

16(1) 

1.311 x 10 10 

1.314 x 10l° 

4 

16(d) 

2.754 x 10 7 

2.831 x 10 7 

4 

16(h) 

4.127 x 10 9 

4.127 x 10 9 

4 

16(m) 

1.874 X 10 10 

1.951 x 10 10 

5 

16(e) 

3.314 x 10 7 

3.654 x 10 7 

5 

16 (i) 

6.342 x 10 9 

6.342 x 10 9 

5 

16 (n) 

2.019 X 10l° 
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TABLE VHI.- NATURAL FREQUENCIES OF A CLAMPED- FREE 120° CONICAL 
FRUSTUM SHELL AS A FUNCTION OF NUMBER OF CIRCULAR 
NODE LINES IN NORMAL DISPLACEMENT w 

[»-*] 


Number of circular 

Mode- identification 

u > 2 , sec -2 

node lines in w 

figure 

Present analysis 

Reference 1 

1 

17(a) 

1.199 x 10 7 

1.195 x 10 7 

1 

17(f) 

1.794 x 10 9 

1.795 x 10 9 

1 

17(k) 

4.290 x 10 9 

4.290 X 10 9 

2 

17(b) 

1.948 X 10 7 

1.953 X 10 7 

2 

17(g) 

1.420 X 10 9 

1.420 x 10 9 

3 

17(c) 

2.420 X 10 7 

2.453 X 10 7 

3 

17(h) 

3.069 X 10 9 

3.076 X 10 9 

3 

17(1) 

8.695 x 10 9 

8.702 x 10 9 

4 

17(d) 

2.940 X 10 7 

3.076 X 10 7 

4 

17 (m) 

1.304 x l 0 l° 

1.325 X 10l° 

4 

17 (i) 

5.224 x 10 10 


5 

17(e) 

3.545 X 10 7 


5 

17 (j) 

1.025 X 1010 


5 

17(n) 

1.760 X 10l° 
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TABLE IX.- NATURAL FREQUENCIES OF A FREELY SUPPORTED SHELL HAVING 
POSITIVE GAUSSIAN CURVATURE AS A FUNCTION OF THE NODAL CIRCLES 

IN NORMAL DISPLACEMENT w 

[»-*] 


Number of circular 
node lines in w 

Mode- identification 
figure 

sec“2 

(present analysis) 

0 

20(a) 

0.1425 

0 

20(f) 

.2340 

0 

20 (k) 

5.659 

1 

20(b) 

0.3607 

1 

20(g) 

3.918 

1 

20(1) 

9.754 

2 

20(c) 

0.5507 

2 

20(h) 

5.914 

2 

20 (m) 

15.90 

3 

20(d) 

0.6641 

3 

20 (i) 

8.571 

3 

20 (n) 

24.40 

4 

20(e) 

0.7260 

4 

20(j) 

12.26 

4 

20 (o) 

35.24 
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TABLE X.- NATURAL FREQUENCIES OF A FREELY SUPPORTED SHELL WITH 
NEGATIVE GAUSSIAN CURVATURE AS A FUNCTION OF THE NUMBER OF 


CIRCULAR NODE LINES IN NORMAL DISPLACEMENT w 

[»-*] 


Number of circular 
node lines in w 

Mode- identification 
figure 

fi2 ? sec"2 
(present analysis) 

0 

21(a) 

0.0246 

0 

21(f) 

1.888 

0 

2 1 (k) 

5.804 

1 

21(b) 

0.1950 

1 

21(g) 

3.243 

1 

21(1) 

8.660 

2 

21(c) 

0.403 

2 

21(h) 

5.128 

2 

21(m) 

13.96 

3 

21(d) 

0.5547 

3 

21(1) 

1.347 

3 

21(n) 

7.738 

4 

21(e) 

0.6499 

4 

21(j) 

11.13 

4 

21(o) 

31.39 
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Figure L- Geometry of a shell of revolution. 
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Figure 2.- Typical idealization of shell of revolution showing geometrically exact finite elements. 
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Figure 3.- Illustration of synthesis of stiffness and mass matrices. 


r = 5 in. = 12.7 cm 
l/Rj = 0 

1/r 2 = 0.2 in." 1 = 0.07874 cm" 1 
r’ = 0 
= 0 

L = 20 in. = 50.8 cm 


E = 

2.96 X 10 7 lb/in2 - 

2.0408 X 10 7 N/cm2 

M = 

0.29 


P = 

7.33 X IQ" 4 lb-se C 2/in 4 = 0.78335 X 10" 4 

h = 

0.008 in. = 0.02032 

cm 


(a) Isotropic cylinder. 


r 

= 1.925 in. = 4.8895 

cm 

l/Rl 

= 0 


1/ R 2 

ss 0.5195 in”* « 0.2045 cm " ^ 

r' 

= 0 


R i 

= 0 


h 

= 0.065 in. = 0.165 cm 

L 

= 15.53 in. = 39.4462 cm 

ph 

= 0.1211 X IQ" 4 lb-sec2/in 3 = 0.3287 X 10 

C 11 

= 1.25 X 10 6 lb/in. = 

= 2.189 X 10® N/cm 

Cl2 

= 0.187 X 10 6 lb/in. 

= 0.327 X 10 6 N/cm 

C 22 

= 0.742 X 10 6 lb/in. 

= 1.299 X 10 6 N/cm 

C 6 6 

= 0.473 X 10 6 lb/in. 

= 0.828 X 10 6 N/cm 

Dll 

= 0.652 X 10® lb-in. 

= 7.367 X 10 6 N-cm 

d 12 

= 1.767 X 10 6 lb- in. 

= 19.964 X 10® N-cm 

d 22 

= 2.767 X 10 6 lb-in. 

= 31.263 X 10® N-cm 

d 66 

= 9 X 10 6 lb- in. = 101.686 X 10 6 N-cm 


(b) Orthotropic cylinder. 

Figure 4.- Properties of the shells analyzed for sample calculations. 
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r = 3 + s \[z/2 (in.) = 7.62 + s ]f3/2 (cm) 

i/r x = 0 

i/Ro = (in _ l) = (cm - l) 

1 0.1667 + s|/3 0.06562 + sl/3 

r' = \j2/2 
Rj = 0 


E = 1 X 10 7 lb/in 2 = 6.8948 X 10 6 N/cm 2 
H = 0.315 

p = 2.54 X 10 -4 lb-sec 2 /in 4 = 0.2714 X 10 -4 N-sec 2 /cm 4 
h = 0.025 in. = 0.0635 cm 

(c) Isotropic 120° conical frustum. 

r = 3 cos(o.5 - - 1.879 (in.) = 7.62 cos(o.5 - - 4.773 (cm) 

l/Rl = 0.333 in -1 = 0.131 cm'l 


R^ = 0 

L = 3 in. = 7.62 cm 
E = 1 lb/in 2 = 0.68948 N/cm 2 
A l = 0.30 

p = 1 lb-sec 2 /in^ = 0.10687 N-sec 2 / cm ^ 
h = 0.001 in. = 0.00254 cm 




(d) Isotropic shell of positive Gaussian curvature. 
Figure 4.- Continued. 


r = 


i/ r i= 

l/ R 2 = 


1 - 20<{cos[-(1.5 - s)0.5] - ly> (in.) 

2.54 - 50.8{cos[-(3.81 - s)0. 01969] - l} (cm) 
-0.05 in- 1 = -0.01969 cm" 1 
-cos[](1.5 - s)0.05] 

-1 + 20{cos[-(1.5 - s)0.05] - l} 

-cos[-(3.81 - s)0. 01969] 


in 


cm 


-1 


r’ = 


-2.54 + 50.8{cos[-(3.81 - s)0.01969] - l} 
sin[-(1.5 - s)0.05] = sin[-(3.81 - s)0.01969] 


Rj = 0 

L = 3 in. = 7.62 cm 
E = 0.91 lb/in2 = 0.6274 N/cm2 
j u = 0.30 

p = 1 Ib-sec2/in4 = 0.10687 N-sec2/cm4 
h = 0.001 in. = 0.00254 cm 

(e) Isotropic shell of negative Gaussian curvature. 


r = 
1/R X = 
l/R 2 = 

r’ = 

R l = 
L = 

E = 

P = 
P = 
h = 


0.5 + s (in.) = 1.27 + s (cm) 

0 

0 

1 

0 

0.5 in. = 1.27 cm 

10.92 lb/in2 = 7.5291 N/cm2 

0.30 

1 Ib-sec2/in4 = 0.10687 N-sec2/cm4 
1 in. =2.54 cm 

(f) Isotropic annular plate. 

Figure 4.- Concluded. 
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2 2 

(Circular frequency) , cu > sec 



Number of circumferential waves, n 


Figure 5.- Minimum circular frequencies of a cylindrical shell computed by present method and method of reference 15. Freely supported edges. 
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(Circular frequency) 2 , a) 2 , sec 



Number of circumferential waves, n 


Figure 6.- Minimum circular frequencies of an orthotropic cylindrical shell computed by present method and method of reference 16. 

Freely supported edges. 
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(Circular frequency ) 2 , cu 2 , sec" 


io 6 r 



Figure 7.- Minimum circular frequencies of a 120° conical frustum by present method and method of reference 1. 

Free-free edge conditions. 
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Figure 8.- Minimum circular frequencies of a 120° conical frustum by present method and method of reference 1. 

Clamped-free edge conditions. 
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, 30 l 1 1 1 1 1 1 

0 2 4 6 8 10 

Number of circumferential waves, n 


Figure 9.- Minimum nondimensional frequencies of a shell of positive Gaussian curvature as computed by present method and 

method of reference 4. Freely supported edges. 
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Number of circumferential waves, n 


Figure 10.- Minimum nondimensional frequencies of a shell of negative Gaussian curvature as computed by present method and 

method of reference 4. Freely supported edges. 
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(Circular frequency) 2 , 



Figure 11.- Minimum circular frequencies of an annular plate as computed by present method and method of reference 17. Free-free edges. 
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(Circular frequency) 2 , a> 2 , sec 


IO"r 



Number of meridional waves, m 

Figure 12.- Frequencies of an isotropic cylinder as computed by present method and method of reference 15. Freely supported edges; n = 2. 


65 


(Circular frequency) , go , sec 


10" 


O 

O 

o ° " ° ° 

CM 

1 

O 

0) 

CO 

CM 

3 


CM 

>N 

u 

c 

<D 

D 

C T ^ 

21 I0 9 

V4— 

o 

"5 

o 


(C 

o 

00 

/ o Calculated values 

/ Faired curve 

10? 

i i i i i 

0 12 3 4 5 

Number of meridional waves, m 


Figure 13.- Frequencies of an orthotropic cylinder as computed by present method and method of reference 16. 

Freely supported edges; n = 2. 
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(Circular frequency) , a) , sec 


10" 


O 

O 

0 n 

1 o 9 

CVJ 1 w 

1 

<J> 

CD 

CO 

o 

\ 2 2 

quency) , cu 

O 

00 


CD 

*4- 

k. 

o 7 

3 I0 7 

O 

w. 

/ 

(C 

o 

CD 

/ o Calculated values 

/ Faired curve 

I0 5 

- / 

I0 4 

_ 1 i i 1 1 l 

0 12 3 4 5 

Number of circular node lines in w 

Figure 14.- 

Circular frequencies of a free-free 120° conical frustum as computed by present method and method of reference 1. 
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2 2 

(Circular frequency) , cu ,sec 


10" 


O 

O 


CSI 

1 

o 

CD 

CO 

CM 

3 


(Circular frequency) 2 

9o 

o Calculated values 
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Figure 15. 

- Frequencies of a 120° conical frustum shell as computed by present method and method of reference 1. Clamped at 

smaller end; n = 2. 
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(a) One nodal circle in w displacement. 



(b) Two nodal circles in w displacement. 




(d) Four nodal circles in w displacement. 


Figure 16.- Natural mode shapes of a free-free 120° conical frustum corresponding to the frequencies shown in figure 14. 
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(e) Five nodal circles in w displacement. 


(f) One nodal circle in w displacement. 


(g) Three nodal circles in w displacement. 


(h) Four nodal circles in w displacement. 


Figure 16.- Continued. 
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(i) Five nodal circles in w displacement. 



(j) One nodal circle in w displacement. 



Figure 16.- Continued. 
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(m) Four nodal circles in w displacement. Five nodal circles in w displacement. 
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(a) One nodal circle in 


displacement 


(b) Two nodal circles in w displacement 


(c) Three nodal circles in w displacement 


(d) Four nodal circles in w displacement 


Figure 17.- Natural mode shapes of a clamped-free 120° conical frustum corresponding to the frequencies shown in figure 15. 
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Figure 17.- Continued. 
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(i) Four nodal circles in w displacement (j) Five nodal circles in w displacement One nodal circle in w displacement (1) Three nodal circles in w displacement 
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(m) Four nodal circles in w displacement (n) Five nodal circles in w displacement 
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Number of circular node lines in w 

Figure 18.- Circular frequencies of a shell with positive Gaussian curvature with freely supported edges, n = 2. 
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Figure 19.- Circular frequencies of a shell with negative Gaussian curvature having freely supported edges, n - 2. 
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(a) Zero nodal circles in w displacement 


(b) One nodal circle in » displacement ( c ) T wo nodal circles in » displacement 


(d) Three nodal circles in w displacement 


Figure 20.- Natural mode shapes of a shell of positive Gaussian curvature corresponding to the frequencies shown in figure 18. 


CO 


W/WMflX U/UMHX 


00 

o 







(e) Four nodal circles in w displacement 


(f) Zero nodal circles in w displacement 


(g) One nodal circle in w displacement 


(h) Two nodal circles in w displacement 


Figure 20.- Continued. 
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(i) Three nodal circles in w displacement 


Cj) Four nodal circles in « displacement 


(k) Zero nodal circles in w displacement 


(1) One nodal circle in w displacement 


Figure 20.- Continued. 
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(m) Two nodal circles in w displacement 


U/UMRX V/VMRX 




(n) Three nodal circles in w displacement 


(o) Four nodal circles in w displacement 


Figure 20.- Concluded. 
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Figure 21.- Natural mode shapes of a shell having negative Gaussian curvature corresponding to the frequencies shown in figure 19. 
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(e) Four nodal circles in w displacement (f) Zero nodal circles in w displacement (fl) One nodal circle in w displacement (h) Two nodal circles in w displacement 
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(i) Three nodal circles in w displacement (j) Four nodal circles in « displacement (k) Zero nodal circles in « displacement (1) One nodal circle in w displacement 




86 


NASA- Langley, 1969 32 


L-6053 


